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In this work we present a model for the propagation of culture on networks of different topology 
and by considering different underlying dynamics. We extend a previous model proposed by Axelrod 
by letting a majority govern the dynamics of changes. This in turn allows us to define a Lyapunov 
functional for the system. 
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INTRODUCTION 

During the last years it has been possible to witness 
the increasing interest in mathematical models aimed to 
describe and analyze social processes. A fruitful sym- 
biosis took place, establishing collaborations among re- 
searchers from both the social sciences as well as physics 
and mathematics. The process gave birth to an inter- 
esting collection of works dealing with a wide spectra 
of social phenomena thus analyzed through a variety of 
mathematical and physical techniques 

Some important and primordial questions constitute 
the motivation to work on these models. It is impor- 
tant to know how the behavior of unorganized individ- 
uals within a society contributes to produce collective 
social phenomena. Next, it is necessary to know how 
stable this emergent phenomena are. Another important 
aspect is knowing to what level social organization is en- 
coded within the topology of the system, or what is the 
extent of the effect of the social structure on the par- 
ticular characteristics of the evolution of a given social 
scenario. A preponderant role regarding this last aspect 
has been played by the social or complex networks. 

Proposing different social networks as schemes for the 
underlying architecture of the society, many authors have 
presented models to describe opinion formation 




rumors diseases |R B and fashion propagation 
urban segregation [lfl| . majority vote etc. In this 
work we will analyze a generalization of a model of cul- 
ture propagation introduced by Axelrod 0,0]. In the 
original model, the cultural background of an individual 
is characterized by a set of F dynamical attributes or 
cultural features, that evolve according to the social en- 
vironment. Each feature, in turn, can take q different 
values, representing possible traits. The individuals are 
located on top of a regular network and interact with 
their neighbors. Through this interaction the cultural 
profile of each individual, and thus the configuration of 
the system, evolve. The interaction is mediated by what 
is called cultural affinity. The more similar an individual 
is to one of its neighbors, the more likely it is that they 
will interact. The interaction consists in the adoption of 
a common trait in one of the F cultural aspects. Typ- 
ically, the system evolves towards a monocultural state, 
but for some parameters values it freezes in a multicul- 



tural state with coexisting spatial domains of different 
cultures. The number of these domains is taken as a 
measure of cultural diversity. 

A systematic analysis of the dependence on q of the 
original model was carried out in [Lj]. Further analysis 
of the model and the role of noise was perform in 
In this work we have modified the original model, al- 
lowing for a wider and more thorough evaluation of the 
cultural environment that envelopes each individual and 
thus influences its cultural tendencies. Each individual 
will evaluate whether or not to copy one of the traits 
adopted by one of its neighbors, basing the decision on 
an observation of the state of its entire neighborhood. A 
sort of majority rule will govern the dynamics of the sys- 
tem. We propose different forms for this majority rule. 
In each function of the state of the system that be- 

haves monotonically in time is found. This function can 
be associated with a Lyapunov function of the system. 
The monocultural state is always the absolute minimum 
of this function, though there exist some local minima, 
corresponding to multicultural situations where the sys- 
tem remains frozen. 



THE MODEL 

Our work will be based on a previous model by Axel- 
rod SE1- We consider that the cultural background of 
any individual can be characterized, in a quite reduction- 
ist way, by a given set of non overlapping features. We 
will call F the number of these features 4>k that define the 
culture in a schematic way. In principle, each feature can 
be associated to a different aspect of the culture such as 
spoken language, preferred foods, or musical style, read- 
ings, sports, etc. In turn, a further subclassification of 
each feature into categories will serve to denote the dif- 
ferent preferences or traits of the individuals regarding 
a particular aspect of his/her cultural background. The 
simplest consideration to achieve such subdivision is to 
consider that any cultural feature may take on any of q 
different values, the same for any <f>k- In principle the 
values are only labels, so despite the fact that we will 
use numbers for classification, any set of symbols would 
work as well. Thus, an individual i will be culturally 
characterized by a cultural vector of F components <f> l k , 
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each one adopting values ranging between 1 and q. The 
way to culturally compare two individuals i and j is to 
measure their cultural overlap, Vu,j) as follows, 

F 

=E<Wi- 
fe=i 

The individuals will be situated on the vertices of a 
graph. Two of them will be considered to be neighbors if 
they are linked by a edge. In the original model, the un- 
derlying network was a bidimensional regular lattice [l2T |. 
The topology of the underlying networks was later gen- 
eralized in [1J|, considering amongst others, Small World 
Networks (SW) 0. In this work, we will consider SW 
networks as well. 

Considering that the individuals will interact with the 
set of their neighbors, it will be useful to define at this 
point the local overlap of a given individual, fli as 

0; = E 

where <r, is the set of neighbors of i, We will also define 
the quadratic local overlap as 

«. 2 = Eir 

Starting from an initial distribution of cultural vec- 
tors, the individuals evolve by analyzing and interacting 
with their environment, adapting their cultural prefer- 
ences according to the tendencies of the neighborhood. 
The original numerical simulations proceed as follows in 
[13, EE EH- At time step t, a randomly chosen individual 
i and one of its neighbors j are evaluated. Their cultural 
overlap is calculated to decide whether they will inter- 
act or not. The interaction takes place with a probabil- 
ity W(Vj)(t) / 'F ', in which case one of the features (f>\ such 
that 4> % k ^ is set equal to 4> J k . Though it is evident that 
W ^ ~ l)i the interaction may affect as well 

the overlaps between i and the rest of the neighborhood 
and thus the change on fij can not be anticipated. 

Interesting results were obtained in E3l by consider- 
ing the whole process of cultural dissemination as an op- 
timization problem. In that work, by analyzing a one 
dimensional system with interaction amongst the first 
neighbors, the authors found a Lyapunov potential that 
allowed them to analyze the stability of the states at 
which the system remained frozen after some evolution 
time. We are interested in extending those results to 
more general situations, namely SWN and other complex 
networks. The global overlap 

i 

can not be claimed to have a monotonic behavior when 
the system, as defined, evolves in time. Suppose that 



as a result of the interaction between i and j at time t 
there is a change in the value of 4>\ . We will call cr™ 1 the 
neighborhood of i such that for any of the fJb individuals 
h rn e cr™ , 4> k m (t) = (j>\ (t) ; and a™ the set of v neighbors 
h n such that <^ n (t) = <j)%{t + 1). The rest of the of 
neighbors will be included in the set a\. If at each time 
step only one change is allowed, we can calculate AS1 = 
fi(t+l)-J2(t) = 2(n i (t + l)-fl l (t)) by considering that 

n 4 (t + i) = X><«)(*)) + (1) 
E - *) + E + !) 

jeo-f jet?? 

= E few)- E i+E 1 - 

jea-i jeo-™- jeer™ 

Thus 

fi i (t+l)=fl f (t)+i/-/i. (2) 

The change in £1 is thus Aft = v — fi, which is not nec- 
essarily equal or greater than zero. By introducing a 
modification in the dynamics of the original model we 
can assure that this condition will be fulfilled and thus 
can think of a Lyapunov functional for the system. 

We will consider different types of dynamics, each one 
associated with a corresponding Lyapunov function but 
at the same time with a clear social interpretation of the 
behavior of the individuals. 

The underlying network will be built up following the 
procedure described in Efl- I n the original model of 
SW networks, a single parameter p, running from to 
1, characterizes the degree of disorder of the network, re- 
spectively ranging from a regular lattice to a completely 
random graph. The construction of these networks starts 
from a regular, one-dimensional, periodic lattice of N el- 
ements linked to 2K nearest neighbors. Then each of the 
sites is visited, rewiring K of its links with probability p. 
Values of p within the interval [0, 1] produce a continuous 
spectrum of small world networks. 

CULTURAL EXCHANGE DYNAMICS 
Case 1: Restricted cultural affinity 

The original model proposed by Axelrod considered a 
very special case of biased dynamics for the interactions 
of the individuals. Despite the fact that individuals are 
immersed in their neighborhood, this was ignored by re- 
quiring that the individual interact with only one of its 
neighbors. Taking this fact into account, a first adapta- 
tion of the original model consists in deciding whether to 
change or not the value of the chosen feature by weight- 
ing the decision with a further evaluation of the influence 
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of the neighborhood. Given that the individual i inter- 
acts with j, the possibility of adopting <f? k for <j> l k will 
depend on the result of an evaluation following a sort of 
majority rule. If by accepting the change of the value 
of <j) l k i will share the value taken by (jf k with a bigger 
group than if by rejecting the change, then i accepts the 
change. With a probability 1/2 the change is accepted 
in case of equality. This is translated into the following 
situation. The change is accepted whenever v > n, and 
with probability 1/2 when v = // . Under this condition, 
Q(t + 1) - Q.(t) > 0. So, we will take £ 1 {t) = -Q,(t) as 
the Lyapunov function of this dynamics, and the system 
will evolve to reach a local or absolute minimum. 



ering sums over cr™, cr™, and a\. On one side we have 

^(t+i) = Efe)W) 2 + ( 3 ) 

^(c (4j) (t)-l) 2 + + 1) 2 

jeer™ jeer™ 

= Efe)W) 2 + E !- 2w w)(*)+ E 1 +Mij)( t ) 

jec, jeer™ 1 jeer™ 

Expanding and regrouping terms we get 

tii(t+i) = ^(t)+^(f+i)-j]r(t+i)+^(t)-ji™(t), 



Case 2: Complete Cultural affinity 

The former rule assures that the global cultural back- 
ground grows or at least is maintained constant, while the 
individuals knows that accepting the change will warrant 
being in bigger group regarding the changing feature. 
But basing the decision of the individual on the com- 
parison of only one feature seems quite myopic. We can 
propose another condition, making the individual base 
the decision on a further evaluation of the local partial 
overlaps, ft™ = £ je<7 m and n ? = Ejea™ w (»,j)- 

Despite that one feature was chosen to be changed, the 
individual decides whether to adopt the new value or not 
by weighting the whole cultural overlap with its neigh- 
borhood and not by analyzing what happens with the 
specific feature to be changed. This is equivalent to say- 
ing that the individual will favor a majority weighted 
by deeper cultural affinity. Now, we can no longer say 
that ft(t + 1) — ft(t) > 0. We must look for another 
quantity £2(t), such that £ 2 (t + 1) — £ 2 {t) < 0. In 
what follows we show that £(t)s = — (ft 2 (i) +ft(i)), with 
ft 2 (i) — j Ej (*)> satisfies the required condition. 

Let us consider that in the proposed interaction be- 
tween i and j, 4> k will be change by <jP k . There are three 
classes of individuals among the neighbors of i, those be- 
longing to cr™, those belonging to <t™ and the rest, that 
will be grouped in a\. The local partial overlap ft- = 
J2jecr' w (i,j) wm 110 De affected, regardless of whether or 
not the interaction takes place. On the contrary, if the 
interaction occurs at time t, ft™(t + 1) = ft™(t) — H and 

n?(t + i) = n?(t) + i/. 

If no interaction is allowed, £<2,(t + 1) = £i(t). If on 
the contrary, the change is accepted, we have 

£ 2 {t + l)-£ 2 (t) = [ft?(t+l)-n?(t)] + [ni(t+l)-fi<(t)]. 



= n?(t) + 2(n?(t + i)-nr(*)) + /*-" 

On the other hand we have Eq. 2. Finally 

Mt + i) - £ 2 (t) = -mut + 1) - nr(*)))- 

The condition to be fulfilled is 

n?(t + i)-fin*)>o 

that corresponds to the imposed constraint. 

It is important to note that in all the cases, the mono- 
cultural state corresponds to the case when the Lyapunov 
function takes an absolute minimum value, £f* . We can 
use this value for normalization, such that Li = £i/£f I . 
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FIG. 1: Proportion of active links p a and of complete overlap 
links pF vs. time, with q = 2 (full), q = 5 (dashed), q = 10 
(dotted), q = 15 (dotted-dashed). Each plot correspond to 
a different value of p: A)p = 0, B)p = 0.01, C)p = 0.5, 
D)p = 0.9. Axelrod' s Case. 



We can expand the rhs of the former equation by consid- 
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NUMERICAL RESULTS 

In what follows we will include results corresponding 
to the cases 1 and 2 as well as those corresponding to the 
Axelrod' s original model, for which we have not defined 
a Lyapunov function. 

We have performed extensive numerical simulations of 
the described model, considering different dynamics. The 
networks have N = 10 4 vertices and connectivity K = 2. 
A typical realization starts with the generation of the 
random network and the initialization of the state of the 
elements. After a transient period, the duration of which 
depends on the parameters of the particular simulation, 
a macroscopic stationary state is achieved. The com- 
putations are then repeated for several thousand time 
steps to perform statistical averages. We consider that 
the system has achieved a macroscopic stationary state 
when the corresponding Lyapunov function of the sys- 
tems reaches a stationary value. We will see that this 
does not imply that the system is steady in a particu- 
lar microscopical state. Indeed, the configuration of the 
system fluctuates among states with equal value of the 
Lyapunov function. There are several aspects character- 
izing the evolution of the system to a stationary value of 
the Lyapunov function. We also analyze the behavior of 
the system subject to the original dynamics, when the 
steady state is not characterized by a Lyapunov function 
but by a microscopically frozen state. 

In the calculations whose results are described in what 
follows, we took F = 10 and several values of q, ranging 
from 2 to 15. At each time step only one change was 
proposed so the system was updated asynchronically. We 
considered that one unity of time corresponded to N time 
steps. 

For each of the dynamics described above, we have 
analyzed several aspects of the evolution of the system. 
First we have calculated the proportion of overlaps p 
between individuals corresponding to three cases, a) po 
when = 0, b)pF when = F and c)p a when 

< 0Juj\ < F. The cases a) and b) correspond to the 
situation when no change in the system is possible be- 
cause the interaction of two individuals: in the case a) 
because no interaction will occur when the cultural over- 
lap is equal to zero, in the case b) because individuals 
are already culturally identical. The only active links 
are those corresponding to he case c). Then we have 
calculated the corresponding Lyapunov function (when 
defined) to show how its value evolves monotonically to 
a steady one. Though this does not provide any infor- 
mation about the inner structure of the system, or about 
the existence of clusters, it helps us to have an idea of 
the amount of cultural differentiation that is present. To 
show that though the Lyapunov function reaches a steady 
value, but the system is not in a steady state, we calcu- 
lated the amount of changes that occur in each time step. 



This was also useful to show that in the Axelrod case, the 
system attained a frozen state, with no changes. 



Axelrod's Case 

This case corresponds to the original model (l2lll4Lll5| 
where individuals interact with only one of their neigh- 
bors at each time step. The interaction is mediated by 
the cultural affinity, defined through the cultural overlap 
■ The stronger the affinity is, the greater the possi- 
bility of interaction between two subjects. In the present 
work, the individuals are located on networks with dif- 
ferent degrees of disorder. The ordered case p — 0, cor- 
responds to a one-dimensional lattice with interactions 
between the first and second neighbors. As stated be- 
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FIG. 2: Proportion of changes vs. time, for different values 
of q and p as in Figure 1. Axelrod's Case. 

fore, we did not find a Lyapunov function for this case, 
and we restrict the displayed results to the time depen- 
dence of the proportion of overlaps pf and p a , and of the 
proportion of changes in the individuals' cultural profiles. 
We recall that po = 1 — Pf + Pa- Figure H displays the 
time evolution, averaged over 1000 realizations, of the 
values pf and p a , corresponding to the amount of over- 
laps ojuj\ — F and < (jJuj) < F, normalized to the to- 
tal number of links, KN. We evaluate these quantities on 
networks with different degree of disorder, namely, p = 0, 
0.01, 0.5,0.9. A deeper insight into what is happening is 
obtained by analyzing the data contained in Figure 
There we show the proportion of changes in the cultural 
vectors of the individuals. Each change corresponds to a 
component of any cultural vector that changed its value. 
We show the amount of changes in a unit of time nor- 
malized to the maximum value allowed N, the number 
of proposed changes. 
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FIG. 3: Proportion of active links p a and of complete overlap 
links pf vs. time, for different values of q and p as in Figure 
1. Case 1. 



FIG. 4: Normalized Lyapunov function L\ vs. time, for dif- 
ferent values of q and p as in Figure 1. 



Figures l.A and 2. A correspond to an ordered under- 
lying network. The system goes to a state where only 
non active links survive, that is, p a — > 0. At the same 
time, while the system reaches a steady state, associated 
with the number of changes approaching 0, the system 
achieves a monocultural state when q < F but goes to 
a multicultural state for higher values of q. These facts 
have been already observed in (l2l Q] . 

When some disorder is introduced into the network, 
the behavior of the system is more complex. By looking 
at FigO we can see that there are two different behav- 
iors for ordered and very disordered networks while the 
intermediate case, p = 0.01 shows a mixture of both. 
The system, in ordered networks evolves rather fast to a 
state of low multiculturality or monoculturality. When 
the disorder is increased, the initial disorder survives for 
longer times. At the end, the system ends in a mono- 
cultural state with except when p = and q > F. Fig- 
ure m showing the number of changes in time confirms 
what was mentioned before. We have not observed sharp 
transitions while varying the p value, the behavior of the 
system indeed, undergoes a smooth change as the spatial 
disorder is increased. 



of the normalized Lyapunov function corresponding to 
this case, L\ is plotted in Figure Q] Starting from the 
ordered case, Figure 3. A, we observe that results do not 
differ so much from what we have previously seen. Again, 
the system goes to a state where only non active links 
survive, reaching a steady state, and achieving monocul- 
turality when q < F and a certain degree of multicul- 
turality for higher values of q. This time we can recur 



10 10 10 10 10 10 10 10 10 10 




Case 1: Restricted cultural affinity 

In the following cases the calculation of the Lyapunov 
function will provide us additional information about the 
system behavior. As in the previous case, Figure dis- 
plays the time evolution, averaged over 1000 realizations, 
of the values pp and p a , evaluating these quantities on 
networks of varying disorder. Figure shows the evolu- 
tion in time of the proportion of changes. The evolution 



FIG. 5: Proportion of changes vs. time, for different values 
of q and p as in Figure 1. Case 1. 

to the Lyapunov function to see that the absolute min- 
imum is reached when q < F, but the system remains 
in a frozen state of multiculturality when q > F. It is 
interesting to observe that the number of changes, Fig- 
ure 03 A, goes to zero. When disorder is included on the 
network, the behavior of the system displays non trivial 
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FIG. 7: Proportion of active links p a and of complete overlap 
FIG. 6: Steady value of L\ vs. p, with q = 2 (full), q — 5 links pf vs. time, for different values of q and p as in Figure 
(dashed), q = 10 (dotted), q = 15 (dotted-dashed). l. Case 3. 



effects as can be can observed in figures 3.B, 3.C and 
3.D. The number of active links is different from zero, 
even when a a steady value for the Lyapunov function is 
reached. Though the monocultural state is the absolute 
minimum, it is note attained by the system, who finishes 
trapped in local minimum. In Figures 0]B, .C and .D 
we observe that the Lyapunov function decreases mono- 
tonically to attain a steady state but not to the absolute 
minimum. On the other hand, the steady values depend 
non monotonically on the disorder of the network. De- 
spite the fact that £\ remains steady, the state of the 
system is not frozen. This affirmation comes from the 
observation of Figure where we find that the number 
of changes remains above zero in all cases. Again, the 
mean value of changes behaves in a non trivial way when 
q or p change. 

Perhaps the most interesting feature is the interplay 
between the effect of the spatial disorder and the values 
of q. This can be better observed by analyzing the be- 
havior of the Lyapunov function. In some cases the disor- 
der introduced by the network prevents the system from 
achieving the previously reached monocultural state, but 
on the other hand, the final degree of multiculturality de- 
pends in a very interesting way from both parameters. A 
new concept should be used now, the concept of plasticity 
of the system. We can link this concept to the amount 
of changes that occur. The plasticity grows with p but 
then decays again. A similar non monotonic behavior in 
terms of p can be observed for L±, plotted in Figure El 



Case 2: Complete Cultural affinity 

The first aspect that we can observe for this case is 
that independently of the degree of disorder of the net- 
work, the state of monoculturality is never achieved, as 
shown in Figure We can again verify the interplay 
between the parameters q and p and their effect on the 
behavior of the system. Another issue to be observed is 
the time scale. The evolution towards a steady value of 
the Lyapunov function is much faster than before. At 
the same time we observe that by increasing the disorder 
the amount of active links grows. This is associated to 
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FIG. 8: Normalized Lyapunov function L2 vs. time, for dif- 
ferent values of q and p as in Figure 1. 

the fact observed in Figure where we can see how the 
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amount of changes in the final state also increases with 
p. The behavior of the Lyapunov function simply veri- 
fies that the system reaches a steady value and that this 
value is far from being the absolute minimum. In all the 
cases, the steady value decreases with q. 




Time 



FIG. 9: Proportion of changes vs. time, for different values 
of q and p as in Figure 1. Case 2. 



CONCLUSIONS 

Axelrod's model shows how a microscopical local pro- 
cess of interaction, leading to convergence provokes the 
emergence of global polarization. In previous works, the 
model was used to analyze the effect of the number of 
cultural aspects and traits on the steady configuration 
of the system. Further analysis of the relative size 
of the largest cultural domain revealed an order disorder 
transition with q, the number of different traits, play- 
ing the role of the control parameter. Under a threshold 
value q c (F), the system converges to a monocultural uni- 
form state. Above q c (F) the system freezes in a multicul- 
tural state, that can be associated to polarization. The 
stability of the multicultural states was analyzed in 0] 
by perturbing the system when frozen in a multicultural 
state and showing the further convergence to the mono- 
cultural sate. Perturbations were associated to cultural 
drift. 

In this work we proposed a different sort of gener- 
alization of Axelrod' s model. We modified the model 
to include interactions among several individuals within 
a neighborhood or to let each individual evaluate the 
changes in its cultural preferences by analyzing those of 
its neighbors. Different ways of considering this extended 
interaction were shown. For each, an associated Lya- 
punov function was found, letting us analyze the conver- 
gence of the system towards an absolute or local minima. 



The disorder of the system was not reduced to that in- 
troduced by the initial condition by increasing the value 
of q, but also included in the spatial distribution of the 
agents. For this purpose we analyzed the effect of the dis- 
order of the underlying network considering small world 
networks of varying disorder. 

The results should be analyzed from several different 
points of view. As expressed in the introduction, we not 
only include initial cultural disorder through higher val- 
ues of q, but also spatial disorder within the underlying 
network. The results linked to this aspects can be com- 
pared with previous results and thus unveil the effect of 
the newly defined interaction of each individual with the 
whole neighborhood. As already known, increasing the 
value of q leads the system to undergo a transition from 
monoculturality to multiculturality. However, the effect 
of spatial disorder attempts against this effect. This can 
be explain by recalling that in a disordered network the 
clusterization of the system is lower and thus, the exis- 
tence of clusters of culture reflected in a polarized situa- 
tion is no longer achieved. When the slightest disorder is 
added to the network, the number of links with overlap 
equal to zero goes to zero. 

The change of the rules of interaction introduces a new 
interesting behavior. Not only does the amount of active 
links not go to zero, with the exception when the under- 
lying lattice is ordered and q — 2, but also the system 
reaches a situation when the Lyapunov function adopts 
a steady value but the system is not frozen. The configu- 
ration of the system changes in time, as can be observed 
from the figures displaying the number of changes in time. 

The results presented here complement what was al- 
ready found in the analysis of the model first proposed by 
Axelrod. The interesting feature is that the system, de- 
spite reaching a steady situation, does not remains static. 
Some aspects still deserve further analysis. Among them 
we will consider in a future work the inclusion of noise. 
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